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Tidal Constraints on Planetary Habitability 



Rory Barnes 1 ' 2 , Brian Jackson 3 ' 4 , Richard Greenberg 5 , Sean N. 
Raymond 2,6 , and Rene Heller 7 

Abstract. We review how tides may impact the habitability of terrestrial- 
like planets. If such planets form around low-mass stars, then planets in the 
circumstellar habitable zone will be close enough to their host stars to experience 
strong tidal forces. We discuss 1) decay of semi-major axis, 2) circularization of 
eccentric orbits, 3) evolution toward zero obliquity, 4) fixed rotation rates (not 
necessarily synchronous), and 5) internal heating. We briefly describe these 
effects using the example of a 0.25 M Q star with a 10 M ffi companion. We 
suggest that the concept of a habitable zone should be modified to include the 
effects of tides. 



1 Introduction 

Exoplanet surfaces are probably the best places to look for life beyond the Solar 
System. Remote sensing of these bodies is still in its infancy, and the technology 
does not yet exist to measure the properties of terrestrial exoplanet atmospheres 
directly. Indeed, the scale and precision of the engineering required to do so is 
breathtaking. Given these limitations, a reliable model of habitability is essential 
in order to maximize the scientific return of future ground- and space-based 
missions with the capability to remotely detect exoplanet atmospheres. 

Here we review one often misunderstood issue: The effect of tides. If the 
distance between a star and planet is small, <0.1 AU, the shape of the planet 
(and star) can become significantly non-spherical. This asymmetry can change 
the planet's orbital motion from that of spherical planets. Simulating the devia- 
tions from the spherical approximation is difficult and uncertain as observations 
of the Solar System, binary stars and exoplanets do not yet provide enough infor- 
mation to distinguish between models. Without firm constraints, qualitatively 
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different models of the planetary response to tides exist. The two most promi- 
nent descriptions are the "constant-phase-lag" and "constant-time-lag" models 
(Greenberg 2009). In the former, the tidal bulge is assumed to lag the perturber 
by a fixed angle, but in the latter it lags by a fixed time interval. Depending on 
the mathematical extension in terms of e, the two models may diverge signif- 
icantly when e>0.3. Throughout this review the reader should remember that 
the presented magnitudes of tidal effects are model-dependent. For more on 
these differences and the details of tidal models, the reader is referred to recent 
reviews by Ferraz-Mello et al. (2008) and Heller et al. (2009). 

We consider tidal effects in the habitable zone (HZ) model proposed by 
Barnes et al. (2008) which utilizes the 50% cloud cover HZ of Selsis et al. (2007), 
but assumes that the orbit averaged flux determines surface temperature (Williams 
and Pollard 2002). We use the example of a 10 M ffi planet orbiting a 0.25 M Q 
star. This choice is arbitrary, but we note that large terrestrial planets orbit- 
ing small stars will be preferentially discovered by current detection techniques. 
This chapter is organized as follows: In § 2 we discuss orbital evolution, in § 3 we 
describe rotation rates, in § 4 we consider the obliquity, and in § 5 we examine 
tidal heating. 



2 Orbital Evolution 

Orbital evolution due to tides should be considered for any potentially habit- 
able world. The asymmetry of the tidal bulge leads to torques which transfer 
angular momentum between rotation and orbits, and the constant flexing of 
the planet's figure between pericenter and apocenter dissipates energy inside 
the planet. These two effects act to circularize and shrink most orbits. In the 
constant-phase-lag model, the orbits of close-in exoplanets evolve in the following 
way (Goldreich and Soter 1966; see also Jackson et al. 2009): 



da _ (G3V GM * R p 2 V^jGjMlRlrn p 



dt V 2 m p Q' p 2 Q' + 



1+^])^ (1) 



de _ ^63 VGMfR 5 p 225 ^/GjMJ%^ \ _ 13/2 



dt V 4 m p Q' p 16 Q'„ 

where a is semi-major axis, G is Newton's gravitational constant, m p is the mass 
of the planet, Q' p is the planet's tidal dissipation function divided by two-thirds 
its Love number, Q'^ is the star's tidal dissipation function divided by two- 
thirds its Love number, R p is the planet's radius, and R* is the stellar radius. 
The Q' values represent the body's response to tidal processes and combines 
a myriad of internal properties, such as density, equation of state, etc. It is a 
difficult quantity to measure, so here we use the standard values of Q'^ = 10 6 
and Q' p = 500 (Mathieu 1994; Mardling and Lin 2002; Jackson et al. 2008a). 
The first terms in Eqs. {H - [2D represent the effects of the tide raised on the 
planet, the second the tide raised on the star. 

Eqs. (Q] - [2]) predict a and e decay with time. As tides slowly change a 
planet's orbit, the planet may move out (through the inner edge) of the habitable 
zone (HZ). This possibility was considered in Barnes et al. (2008), who showed 
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Figure 1. Contours of equilibrium rotation period in days for a 10 M ffi 
planet orbiting a 0.25 M star. The gray region is the HZ from Barnes et 
al. (2008). 

for some example cases the time for a planet to pass through the inner edge of 
the HZ. Such sterilizing evolution is most likely to occur for planets with initially 
large eccentricity near the inner edge of the HZ of low mass stars (<0.3Mq). 
Even if a planet does not leave the HZ, the circularization of its orbit can require 
billions of years, potentially affecting the climatic evolution of the planet. 

3 Rotation Rates 

Planetary rotation rates may be modified by tidal interactions. Although planets 
may form with a wide range of rotation rates f2, tidal forces may fix O such that 
no net exchange of rotational and orbital angular momenta occurs during one 
orbital period. The planet is then said to be "tidally locked," and the rotation 
rate is "pseudo-synchronous" or in equilibrium. The equilibrium rotation rate 
in the constant-phase-lag model is 

19 

n eq = n(l + -e 2 ), (3) 

where n is the mean motion (Goldreich 1966). Note that planets only rotate 
synchronously (one side always facing the star) if e = (the constant-time-lag 
model makes the same prediction). Therefore, the threat to habitability may 
have been overstated in the past, as independently pointed out by several recent 
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Figure 2. Time in years for a 10 M© planet orbiting a 0.25 M© star to 
evolve from an obliquity ip — 23% to 5? The gray region is the HZ from 
Barnes et al. (2008). 



investigations (Barnes et al. 2008; Ferraz-Mello et al. 2008; Correia et al. 2008). 
Figure [1] shows the values of the equilibrium rotation period for a 10 planet 
orbiting a 0.25 M star as a function of a and e. 



4 Obliquity 

Tidal effects tend to drive obliquities to zero or tt. The constant-time-lag model 
of Levrard et al. (2007) found a planet's obliquity ip changes as 

chp = sin (VQifp / cos (Vpe^ _ g \ 
dt Cpfio^ \ n / 

where 



and 



1 + 3e 2 + |e 4 

Cl ~ (1 _ e 2)9/2 ' ^ 
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Figure 3. Tidal heating fluxes for a 10 M® planet orbiting a 0.25 Mq star. 
Contour labels are in W m~ 2 . The dashed contours represent the boundaries 
of the tidal habitable zone (Jackson et al. 2008c; Barnes et al. 2009b). The 
gray region is the HZ from Barnes et al. (2008). 

In the preceding equations r g)P (= 0.5) is the planet's radius of gyration (a 
measure of the distribution of matter inside a body), Qq is the initial rotation 
frequency, and r p is the "tidal time lag", which in this constant-time-lag model 
replaces Q' p . We assumed Q' p = 500 for the planet at its initial orbital configura- 
tion and set r p = l/(nQ' p ), i.e. initially the planet responds in the same way as 
in a constant-phase-lag model. In the course of the orbital evolution, r p was then 
fixed while n and Q p evolved in a self-consistent system of coupled differential 
equations. In Fig. [2] we show the time for a planet with an initial obliquity of 
23?5 to reach 5°, a value which may preclude habitability (F. Selsis, personal 
communication). However, obliquities may easily be modified by other planets 
in the system (Atobe et al. 2004; Atobe and Ida 2007) or a satellite (Laskar et 
al. 1993). 



5 Tidal Heating 

As a body on an eccentric orbit is continually reshaped due to the varying 
gravitational field, friction heats the interior. This "tidal heating" is quantified 
in the constant-phase-lag model as 



H = 



63 (GM«) 3 / 2 M«itf 
"4" % 



(9) 
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(Peale et al. 1979; Jackson 2008b). However, in order to assess the surface effects 
of tidal heating on a potential biosphere, it is customary to consider the heating 
flux, h = H/4irRp, through the planetary surface. Jackson et al. (2008c; see also 
Barnes et al. 2009b) argued that when h > 2 W mT 2 (the value for Io [McEwen 
et al. 2004]) or h < 0.04 W m~ 2 (the limit for plate tectonics [Williams et 
al. 1997]), habitability is less likely. Barnes et al. (2009b) suggested that these 
limits represent a "tidal habitable zone" . In Fig. [3] contours of tidal heating are 
shown for a 10 planet orbiting a 0.25 Mq star. The tidal habitable zone is 
the region between the dashed curves. Note that a and e evolve as prescribed 
by Eqs. (Q]-[2]), and hence the heating fluxes evolve with time as well. 
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